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Abstract 

It is possible to associate temperatures with the non-extremal horizons of a large class of 
spherically symmetric spacetimes using periodicity in the Euclidean sector and this procedure 
works for the de Sitter spacetime as well. But, unlike e.g., the black hole spacetimes, the de 
Sitter spacetime also allows a description in Friedmann coordinates. This raises the question 
of whether the thermality of the de Sitter horizon can be obtained, working entirely in the 
Friedmann coordinates, without reference to the static coordinates or using the symmetries of 
de Sitter spacetime. We discuss several aspects of this issue for de Sitter and approximately 
de Sitter spacetimes, in the Friedmann coordinates (with a time-dependent background and the 
associated ambiguities in defining the vacuum states). The different choices for the vacuum states, 
behaviour of the mode functions and the detector response are studied in both (1+1) and (1+3) 
dimensions. We compare and contrast the differences brought about by the different choices. In 
the last part of the paper, we also describe a general procedure for studying quantum field theory 
in spacetimes which are approximately de Sitter and, as an example, derive the corrections to 
thermal spectrum due to the presence of pressure-free matter. 

1 Introduction 

There exists an extensive literature on quantum field theory in de Sitter spacetimes (for a non 
exhaustive sample, see Ref. [HE])- From a theoretical point of view, the high level of symmetry 
exhibited by the de Sitter geometry makes it an important and tractable example. On the other 
hand, observations suggest that the evolution of our universe is described by a (near) de Sitter 
geometry both during the early inflationary phase as well as during the current accelerated phase 
of expansion. While quantum eff'ects are not expected to play a serious role in the current phase 
of the expansion (see, however, [3]), they play an important role during the inflationary phase and 
possibly seeds the cosmic structure we see today. This was part of the motivation to study quantum 
field theory in de Sitter background. 
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A key feature which arises in such a study of the de Sitter spacetime is the thermal nature of the 
vacuum state and the temperature which one can associate with the horizon. One of the purposes of 
this paper is to explore aspects of this thermality from different perspectives. Since it is a fairly well 
trodden path, we will begin by describing the motivation (for yet another paper!) and the specific 
point of view adopted here. 

The connection between spacetime horizons and thermodynamics is quite well-known. The first 
example of this kind, of course, was in the case of black hole spacetimes [1]. This was soon followed by 
the discovery of horizon temperature in much wider contexts like, for example, in the case of Rindler 
spacetime [5] and de Sitter spacetimes [B]. While all these three spacetimes (black hole, Rindler, de 
Sitter) have very similar metric structure — when one uses static coordinates which cover part of the 
manifold — they also have significant differences [7]. For example, one can provide a fairly rigorous 
geometrical description of the black hole horizon. In contrast, the horizon in the Rindler spacetime 
is just a null surface in flat spacetime as viewed by uniformly accelerated observers and is clearly 
observer dependent. The situation as regard de Sitter spacetime is somewhat in between: While 
it has a geometric description, it also shares with the Rindler spacetime certain level of observer 
dependence. The location of the horizon surface will be different even for two different observers 
who are translated with respect to each other by a spatial vector. More importantly, the de Sitter 
spacetime allows, in addition to the static coordinate system, a spatially homogeneous Friedmann 
coordinate description. In these coordinates, quantum field theory reduces to the quantum mechanics 
of an oscillator with time dependent frequency with the well known difficulties arising from the 
ambiguity of vacuum state in a time dependent background. The black hole spacetime, of course, 
does not posses such a spatially homogeneous description. The Rindler spacetime does but, in this 
case, the homogeneous description is just the inertial coordinates in flat spacetime with no time 
dependence — thereby making that description trivial. 

Another issue which crops up in the study of de Sitter spacetime is the following. Neither 
the inflationary phase nor the currently accelerating phase of the universe is strictly de Sitter. 
While one can provide a fairly elegant mathematical description of quantum field theory in exact 
de Sitter spacetime, many of these techniques will fail when the spacetime is only approximately de 
Sitter. It is interesting to ask how much progress one can make in studying such (approximately 
de Sitter) spacetimes and how much of the results valid in exact de Sitter will continue to hold 
(in an approximate sense) in such spacetimes. For example, it is fairly complicated to work with 
an analogue of approximately static coordinate system when the universe is not strictly de Sitter. 
Many of the techniques used to define the vacuum states in exactly de Sitter spacetimes will also be 
inapplicable when the manifold has no de Sitter symmetry. 

It is thus clear that, given the special features possessed by de Sitter spacetime, one can approach 
the problem of quantum field theory in de Sitter spacetime from many different perspectives, not 
all of which will be easily generalizable to an approximately de Sitter spacetime. This motivates 
us to examine closely several aspects of quantum field theory in de Sitter spacetime delineating the 
properties which arise (in one way or the other) from the symmetry of de Sitter spacetime from those 
which are of more general nature. Such a study also reveals some significant differences between de 
Sitter spacetime in (1+1) dimension and de Sitter spacetime in (1+3) dimensions. 

We will now briefly describe some of these issues which will be discussed in detail in the paper. 
As we mentioned earlier, there is a very standard procedure for obtaining the thermality of the 
horizon in the static coordinate system. This procedure works for a very wide class of spacetimes 
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and can, for example, handle Rindler, black hole and de Sitter spacetimes at one go. But when 
the cosmological spacetime is not exactly de Sitter, no static coordinate system will exist. One can 
still define an "approximately static" coordinate system but this proves to be difficult to handle 
mathematically. 

If one decides to work with a Priedmann coordinate system, then the mathematics simplifies 
considerably because we will be dealing with a quantum mechanical problem rather than quantum 
field theory. But conceptually, we now have to tackle the issue of defining the vacuum state in a 
time dependent background. This turns out to be reasonably straight forward in (1+1) dimension 
in which conformal invariance of a massless scalar field helps the analysis. But in (1+3) dimension 
it is not possible to have sensible limits for the mode functions in the infinite past if one works with 
the massless scalar field $(t,x) as the primary variable. The usual trick is to work instead with the 
variable x(iix) = a(t)$(t,x) and define a vacuum state in the asymptotic past for xit^^)- When 
t — )• — oo, a{t) — )• and it is the scaling out of this factor which allows us to define a sensible vacuum 
state in (1+3) dimensions. The resulting vacuum state is the well known Bunch-Davies vacuum [8] 
for the de Sitter spacetime. Clearly, this depends on the behaviour of a{t) and there is no natural 
analogue of this vacuum state for non-de Sitter spacetimes. 

An alternative to the above procedure is to define a vacuum state at some fixed time t = to, say, 
by choosing the modes which behave as close to the positive frequency modes as possible at this 
instant. We will call this the co-moving vacuum since it is based on the co-moving time coordinate 
of Friedmann spacetime. In general, this vacuum state differs from the Bunch-Davies vacuum but 
it has the advantage that the evolution of a{t) for t < to becomes irrelevant for its definition. It is, 
therefore, well suited to study spacetimes which are de Sitter at late times with deviations from de 
Sitter geometry in the early epochs. 

Once the vacuum state is defined, in the asymptotic past or at some other chosen moment, one 
could study the mixing of positive and negative frequency modes due to the time dependence of the 
background expansion. In particular, one would be interested in knowing whether the mixing leads 
to a thermal nature for the state at later times. It does happen in the case of (1+1) dimension but 
the spectrum is not strictly Planckian in the case of (1+3) dimensions. There are some interesting 
peculiarities which arise in this context when we try to obtain thermality working entirely in the 
Friedmann coordinates. 

Finally, one can also study the inter-relationship between the mode functions defined in Fried- 
mann coordinate system and those defined with static coordinates. This is an exercise in evaluating 
the Bogolioubov coefficients and we do find that one recovers standard Planck spectrum without any 
deviation. This allows us to establish a correspondence between the vacuum states defined using 
the two coordinate systems but — since static coordinate systems do not exist for approximately de 
Sitter spacetimes — the approach does not allow an easy generalization to more realistic cases. 

In the last part of the paper we study the mode functions in approximately de Sitter spacetimes 
in Friedmann coordinates. We find an explicit solution to the wave equation, correct to the necessary 
order of approximation, and use it to describe the deviations from the exact de Sitter spacetime. 
This approach is quite general and is capable of handling a wide variety of cases when the evolution 
is approximately de Sitter. 

The plan of the paper is as follows: We briefly review thermal aspects of horizons in static 
coordinate syatem in Section [2j In Section [3l we solve for the modes of a massless scalar field in 
spatially flat de Sitter spacetime in (1+1) and (1+3) dimensions and define the Bunch-Davies and 
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comoving vacuum states in the Priedmann coordinate patch. These modes evolve in time and the 
physical content of the modes at later times is determined by evaluating the mixing coefficients in 
Section 21 working entirely in the Friedmann coordinates. In Section [5l we study the response of a 
detector coupled with the field as a way to provide an operational meaning to the mixing coefficients. 
We next compare, in Section[6l the mode functions defined in the Friedmann patch with those defined 
in the static patch of the de Sitter spacetime to reproduce some standard results. Finally, we study 
the corresponding effects in the quasi-de Sitter geometry. We consider a small perturbation to the 
de Sitter metric and develop the perturbative framework to find the corrections to the field modes 
and the corresponding power spectrum in the quasi-de Sitter case. This procedure is illustrated by 
taking the model of the universe containing pressure-free matter and cosmological constant which 
behaves like quasi-de Sitter at late times. Section [8] describes the conclusions. 



2 Thermality in static coordinates 

To set the stage, we shall begin by briefly reviewing some well known results (see e.g., chapter 14 
of Ref. [9] for more details) related to the temperature of horizons in static coordinates. Several 
spacetimes of interest including the Schwarzschild, de Sitter, Rindler, can be described by a line 
element of the form 

ds^ = f{r)dt^ --^^-dLl (1) 

where dL^ is the transverse metric and /(r) vanishes at the horizon, r = ro with f'{rQ) = 2k ^ 0. 
Then, using a Taylor series expansion near the horizon, we can write / ~ 2kI where I = {r — tq) and 
the metric near the horizon takes the form: 

ds^ = 2kI dt^ - — - dL\. (2) 

In the case of Rindler spacetime, this is exact and k denotes the acceleration of the Rindler observer. 
In other cases, the metric reduces to this form close to the horizon with k denoting the surface 
gravity. 

This (Rindler) form of the metric makes it obvious that the singular behaviour of the metric near 
/ = is a coordinate artefact. It is possible to introduce several, different, sets of coordinates which 
will cover the entire manifold without any pathology at the horizon. One such choice, (T, X), which 
we will call Kruskal-like coordinates is obtained by the transformations: 

f dr 

= e'"'* coshKt; kT = e'"'* sinh r^= / 77^ (3) 



which lead to the metric 



ds'' = ^,^ {dT^ - dX^) + dLl (4) 

that covers the full manifold. Here / should be treated as a function of {T,X). The horizon at 
r = ro IS now mapped to = but with the factor //(X^ - T^) remaining finite at the horizon. 

It is now possible to show that the vacuum state of a quantum field defined on the T = 
hypersurface appears as a thermal state to observers confined on the right wedge X > \T\. This is 
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most easily done by making an analytic continuation to the imaginary time coordinates by Te = —iT 
and tE = —it. The time evolution of the system in terms of Te will take the field configuration 
from Te = to Te — )• oo and will be governed by a global Hamiltonian Hgi. One can equivalently 
describe the same evolution in terms of t^;, which behaves like an angular coordinate from t^; = 
to tE = 27r/K when we use the Hamiltonian Hgt which determines time evolution in the static time 
coordinates. The entire upper half-plane T > can be covered in two completely different ways: 
either in terms of the evolution in Te or in terms of the evolution in tE- In (Te,X) coordinates, we 
vary X in the range (—00,00) for each Te and vary Te in the range (0, 00). In {tE,x) coordinates, 
X varies in the range (0, 00) for each tE which varies in the range (0, tt/k) like an angular variable. 
This allows us to prove, using standard path integral techniques [9l[T0] that, 

{vacl^L, (Pr) « (0L|e""^-/l0R). (5) 

where (pL and are the field configurations in left and right parts of the plane on the T = 
hypersurface. One can find the density matrix for observations confined to the right wedge by 
tracing out the field configuration (f)^ on the left wedge. This computation gives: 

which is thermal with the temperature f3~^ = K,/2Tr. Thus, the vacuum state of the field defined on 
the T = hypersurface leads to a thermal density matrix with temperature k/2tt as far as static 
observers in the right hand wedge are concerned. 

In the case de Sitter the metrics in the static and Kruskal-like coordinates 

are connected by the coordinate transformations: 

From the form of the metric in in the two coordinate systems, it is obvious that the (T, X) coordinate 
system is not static because the metric depends on T. On the other hand, the coordinate system 
(t, x) which covers the right wedge has a static time coordinate t. It is well known that defining 
a vacuum state in a time dependent background is non-trivial and often ambiguous. In the above 
analysis we have chosen to define a vacuum state at a particular space-like hypersurface T = and 
examine its properties in terms of the static coordinates. A different definition for the vacuum state, 
in general, will lead to a different description in static coordinates. We will see that similar issues 
arise later on when we study de Sitter universe in the Priedmann coordinates as well. 

An alternative procedure to determine the thermal nature of the horizon is based on the cal- 
culation of relevant Bogoliubov coefficients. Since we have two coordinate systems — Kruskal-like 
and static — covering part of the manifold, one can obtain, in principle, the relation between the 
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field modes which are natural to these coordinate systems and compute the Bogoliubov coefficients 
between them. Let the field modes be given by some functions (j){T, X) and "f") in the Kruskal-like 
and static coordinates respectively in the region of the manifold where both are well defined. (For 
simplicity, we have ignored the dependence on the transverse coordinates which play no role in the 
discussion, as we shall see.) It is often not possible to obtain closed expressions for the field modes 
due to mathematical complexity. However, it is possible to evaluate the Bogoliubov coefficients 
using a simple trick: Since the Bogoliubov coefficients that relate the two sets of field modes are 
independent of the hypersurface which is used to evaluate the Klein-Gordon inner product, we can 
choose this hypersurface to be arbitrarily close to the horizon. The field equations reduce to a two- 
dimensional wave equation near the horizon making the dependence in the transverse coordinates 
(and the mass of the field) irrelevant. Conformal invariance then allows us to determine the field 
modes near the future horizon which take the form of plane waves in the relevant coordinates. That 
is, 

X.(t,r) = -^e— (9) 

and 

</,fc(r,X) = -l=e-'=^. (10) 



where u = t — r* and U = T — X respectively. These are related by kU = — e""" which signifies an 
exponential redshift near the horizon. As is well known, the relevant Bogoliubov coefficient (which 
we will have the occasion to evaluate explicitly later on) will now lead to a thermal spectrum of 
particles. 

This discussion is, of course, applicable to the de Sitter universe described by the metric in Eq. ([1]) 
with /(r) = 1 — H^r^ and will lead to a temperature H/2-k. More precisely, if we introduce Kruskal- 
like coordinates in the de Sitter manifold and define a vacuum state on the T = hypersurface, 
then such a vacuum state will lead to a density matrix with temperature H/2tt for the observers 
confined to the region r < H~^. Once again, it should be stressed that the de Sitter metric in the 
Kruskal-like coordinates is not static and the vacuum state is defined using the T = hypersurface. 

This analysis is completely in tune with what could be done in black hole spacetimes as well 
as in the case of Rindler spacetime. But in the case of de Sitter we have an alternative coordinate 
system available to us, viz.., the standard Friedmann coordinate system. This allows us to study the 
dynamics of a quantum field entirely in the Friedmann coordinate system and explore whether we 
can recover the thermality of the horizons and other features. In such a study we necessarily have 
to work with a time dependent background but — as we have emphasized above — this is implicit 
even when we use Kruskal-like coordinates and relate them to static coordinates. We can, therefore, 
adopt a similar strategy in the Friedmann coordinate system by defining a vacuum state at some 
suitable hypersurface and studying its particle content as the evolution proceeds. 

This approach has one extra advantage. The static coordinate system exists only for the exact de 
Sitter universe. When there are deviations from de Sitter nature, we can still describe the universe 
in a very natural fashion using Friedmann coordinate system. But in this case, we will not have the 
luxury of an alternative static coordinate system to describe the physics. Therefore, a formalism 
which addresses issues like thermality working entirely in Friedmann coordinate system, without 
using any of the symmetries of the de Sitter universe, is well suited for the study of near de Sitter 
geometry. We will find that obtaining thermal nature of the horizon working entirely in Friedmann 
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coordinates is — surprisingly — not an easy task. In fact we could not find any previous work 
in published literature which discusses quantum field theory in de Sitter spacetime from such an 
approach and obtains thermal nature of the horizon. We shall now turn to this study. 



3 Massless scalar field modes in de Sitter spacetime 

Throughout the paper, we will confine ourselves to massless, minimally coupled scalar field in de 
Sitter spacetime. The action for the field $(t,x) is given by: 

Sm = \j (11) 

It turns out that the dynamics is somewhat different in (1+1) dimensional spacetime compared to 
(1+3) dimensional spacetime. We will first study the behaviour in (1 + 1) and then follow the same 
procedure for the (1 + 3) case. This will bring out the similarities and some curious differences 
between the two cases. 

3.1 Field modes in dS2 spacetime 

We will describe the de Sitter spacetime in Priedmann coordinates with k = Q. Then, the (1+1) 
dimensional metric is given by, 

ds^ = dt^ - a\t)dx^ (12) 

and the field equation reads, 

d?<P + -dt<^ - \dl'i> = 0. (13) 
a 

We decompose $ in terms of a complete set of orthonormal functions /k in the form 

^{t,x)= / — afcA + 4/,* . (14) 



j-oo2vr . 

Spatial homogeneity allows us to separate out the x dependence and write: 

A(x,t) = e^^^V|fc|(t) (15) 
Substituting in Eq. (fT3|) and solving the resulting equation, we find that 



V'fc(t) = Asfc(i) + B^sl(t); Skit) = exp l-ik [ 4^ ) (16) 

(The k in these expressions actually stand for \k\; we will not explicitly show the modulus sign 
hereafter for notational simplicity.) The result is obvious from the fact that in (1+1) dimension the 
scalar field action is conformally invariant and any Friedmann spacetime is conformally flat with the 
conformal time coordinate rj defined through dr] = dt/a(t). For dS2 with the scale factor a{t) = e^*, 
the solution is, 



..(t) = -l=exp 



(17) 
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where the phase ensures that, in the H ^ hmit, the mode function reduces to standard flat 
spacetime modes. The constants Ak and in Eq. (jl6p are determined using the appropriate 
boundary conditions and thus decide the choice of the vacuum state for the field. For example, when 
i7 — )• 0, the choice Ak = 1 and 8^ = gives the positive frequency mode and selects the standard 
inertial vacuum for the flat spacetime. 

In the presence of an expanding background, it is difficult to define a unique choice for the 
vacuum and we need to study different choices and their physical properties. One possible choice 
would be to define the vacuum state at the asymptotic past by choosing the field modes such that 
they reduce to positive frequency modes in this limit. It is, however, clear that the mode function 
in Eq. ()17p does not have a well defined phase when t — t- — oo. (This is related to the fact that 
in the asymptotic past a — )• 0.) The usual procedure adopted in the literature to circumvent this 
problem is to abandon the idea of defining a vacuum state using the t coordinate and instead use 
the conformal time rj which, for the de Sitter universe, can be taken to be r/ = (I - e~^^)/H. (The 
integration constant is chosen to give the correct limit of r/ — t- t when H ^ 0.) Then our mode 
function in Eq. (jl7p 

is indeed a positive frequency solution with respect to rj (at all times) and therefore the choice Ak = 1 
and Bk = gives a natural choice for the vacuum. This is the conventional Bunch-Davies vacuum 
defined with respect to conformal time by the choice of mode functions 



V'i^'^^ (t) = exp 



2k 



(19) 



While the Bunch-Davies vacuum is the preferred choice in the literature, it is clear that it is 
more in tune with the conformal time coordinate r] than with the co- moving time coordinate t. In 
the Friedmann metric, the co-moving time t has a direct physical significance as the proper time of 
the co-moving, geodesic clocks. This motivates us to look at the possibility of defining a co-moving 
vacuum with mode functions which behave as close as possible to the positive frequency modes with 
respect to co-moving time coordinate t. We can take a cue from the discussion in the last section 
where we saw that, even in the Kruskal-like coordinates for the de Sitter spacetime, the metric is 
time dependent and the vacuum state is defined on a particular hypersurface T = 0. In a similar 
fashion, we can choose the modes in Eq. ()16p by demanding that at some time t = to they behave 
like positive frequency modes. Because of the time translational invariance, we can take to = 0, 
without the loss of generality, as long as to is finite. That is, we impose the conditions: 

MO) = -^e-^%=o ; 4(0) = ^ e-^%=o. (20) 
\' 2k \J 2k 



(The same physics is obtained if we take to 7^ with the replacement of k by ke~^^" which ensures 
that k is the co- moving wave number defined at to = 0). These conditions imply that at t = to(= 0) 
the mode function and its derivative behave like a positive frequency mode. 

We can now determine the coefficients Ak and Bk using this condition and — somewhat curiously 
— we will again find that Ak = 1 and Bk = 0. That is, the mode function ijj'k^^^ {t), evolved from the 
co-moving vacuum defined at t = to(= 0) is same as Bunch-Davies state '4'k^^\'t) defined earlier in 
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dS2. This result is independent of the choice for to thereby showing that the Bunch-Davies vacuum 
can also be interpreted as a co-moving vacuum state defined using the conditions in Eq. (j20p . 

We will see later, this equivalence is a special feature of (1+1) dimension and does not hold in 
(1+3) dimensions where the co-moving and Bunch-Davies vacua are different. 



3.2 Field modes in dS4 spacetime 

We shall now follow the same procedure as above in the (1 + 3) dimensions. The metric is now given 
by 

ds"^ = dt^ - exp{2Ht)dx^ (21) 

where H is the Hubble constant and sets the only length-scale (or time-scale) in the problem to be 
1/H. The field equation for $(t,x) in this metric reads: 

d^^ + 3Hdt^ -exp{-2Ht)dl^ = (22) 

As usual we expand the field in terms of a complete set of orthonormal functions /k and write: 

$(x,t) = I ^{ak/k(t,x) + aU(t,x)} (23) 

where spatial homogeneity allows us to express the field modes in the form: 

fUt,^) = e'''-^Mt) (24) 
where k = |k|. The equation in ipkit) then becomes, 

■ipk + SHipk + exp(-2m)fcVfe = 0. (25) 

with the solution, 

Mt) = ^kSkit) + BksUt) (26) 

where 



V2fc 



f + e"^*) (27) 



(Because of the existence of Ak and Bk, the normalization of Sk is not unique; we choose it in 
such a way that, when Ak = 1 and = 0, the functions satisfy the standard orthonormality 
conditions with respect to the Klein-Gordon inner product.) Again, the constants Ak and Bk are 
to be determined using the appropriate boundary conditions which makes a choice for the vacuum 
state for the field. In (1+3) dimensions also, we see that when — the choice Ak = 1 and Bk = 
leads to the standard positive frequency mode in flat spacetime and selects the inertial vacuum. Our 
interest is to explore the different choices in the presence of expanding background. 

As in dS2 case, let us first study the behaviour of the modes in the asymptotic past. We see 
that, in the t — )• — oo limit, the expression in Eq. ([27|) goes to: 

-Ht\ -Ht 



.,(i)^-^exp^-e-^*)e-^*. (28) 
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This does not have a well-defined limit and hence cannot be used to define a vacuum state for the 
field. In this respect, both (1+1) and (1+3) dimensional results are similar. 

We found that, in the (1+1) dimensional case we could use the conformal time coordinate r] to 
define a natural vacuum state in the asymptotic past. In the present case, however, the situation is 
different. In terms of the conformal time r], the mode function becomes (in the asymptotic past): 

Sk{r]) ^ (29) 

and we see now the crucial difference from the (1 + 1) dimensional case. There is an extra a{r]) in 
this case which prevents us from treating it as the standard positive frequency mode. 

The result also suggests a possible way-out which is usually adopted in the literature. Instead of 
quantising $ we may choose to quantise $ = a(t)$. This is a (time dependent) point transformation 
of the dynamical variable which is permissible in the classical description. We then see that the 
choice Ak = 1 and Bk = will give the modes exp(— i/cry) which has the standard form for ^, when 
treated as a quantum field. This is the usual procedure in the literature and this choice leads to 
the conventional Bunch-Davies vacuum. But note the the situation was not as straightforward as in 
the case of (1 + 1) dimensions and we needed to remove a factor a{t) to define the vacuum state in 
(1+3) dimensions. 

The difference is more acute when we try to define a co- moving vacuum. As in the case of (1 + 1) 
dimension one can define the co-moving vacuum by imposing the conditions given in Eq. (|2U|) and 
thus determining At and Bk- Because of the time translation invariance, we can again define the 
co-moving vacuum at at t = and the result for any other time, to can be obtained by a finite shift. 
Hence the conditions we impose on the modes are 

MO) = -4^^~'%=o ; UO) = ^ e~^%=o. (30) 



These allow us to determine the constants Ak and Bk as: 

H + 2ik ^ H 

which define the mode function, Tp^^'^^\t), evolved from the co-moving vacuum choice defined at 
t = 0. 

When we did this in (1+1) dimension, we found that Ak = 1 and Bk = - instead of the 
expressions in Eq. (|3ip — thereby showing the equivalence of co-moving and Bunch-Davies vacuum. 
But in (1+3) dimensions we get a different result, viz. that the co-moving vacuum is different from 
the Bunch-Davies vacuum. The difference can be traced, algebraically, to the existence of the a(ry) 
factor in Eq. (|29D . 

To summarise, we can define the vacuum states by imposing suitable boundary conditions on the 
mode functions and thus determining the constants Ak and Bk- If we work in the asymptotic past, 
then one can choose the modes to be exp{—ikr]) in (1+1) dimension, thereby defining the Bunch- 
Davies vacuum. In (1+3) dimensions this is not possible with the original scalar field. But if we 
work with a{t)^, instead of one can again define the modes such that they behave as exp{—ik7]) 
in the asymptotic past. Alternatively, one can attempt to define a co-moving vacuum by imposing 
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the condition that the modes must behave as close to positive frequency solutions as possible, with 
respect to the co-moving time coordinate t, at some time t = to. Because of time translation 
invariance, we can choose to = without loss of generality. We then find that, in (1+1) dimension, 
the co-moving vacuum is equivalent to the Bunch-Davies vacuum. But in (1+3) dimensions, these 
two mode functions (and hence the vacua are different). We shall now explore the properties of these 
vacuum states. 



4 Evolution and mixing coefficients at later times 

The Bunch-Davies and the co-moving vacua are defined by the condition that the mode function is 
purely positive frequency at a given moment of time t = to. In the case of Bunch-Davies vacuum, 
this is done in the asymptotic past {to — >• — oo) while in the case of co-moving vacuum we choose 
this to be to = 0. Once this initial condition is set, expansion of the universe will evolve the mode 
functions to a mixture of positive and negative frequency modes, with respect to the co-moving time 
coordinate, at any later time. This mixing can be analysed in terms of two mixing coefficients, ai, 
and in the expansion: 

Mt)= / — (a.e"-* + /3,e-*) (32) 



It is slightly more convenient to let the frequency vary over both positive and negative values and 
write: 

/°° d;/ 
^^/(^)e-^^* (33) 

so that 

a, = f{u), p, = f{-v)- z.>0 (34) 

The mixing coefficients defined by Eq. (j32p are similar to Bogoliubov coefficients but not the same. 
We stress that in Eq. (|32p ipk{t) is expanded in terms of the complete set of orthonormal functions 
exp(±iz/t) which are not the solutions to scalar field wave equation in the de Sitter background. 
Physically, one can think of these functions exp(±iz/t) as defining the instantaneous positive and 
negative frequency mode functions with respect to the co-moving time. But as we shall see, these 
mixing coefficients have interesting properties and in fact play a direct role in the response of detec- 
tors. We shall say more about it later on. 

The task of determining the mixing coefficients is thus reduced to calculating the the Fourier 
transform of ipkit), 

/oo 
dte^'^'Mt)- (35) 
'OO 

Often we will be interested in |a^p and which can be obtained from the power spectrum |/(z^)p. 
We shall now compute these for the different cases. 

4.1 Mixing coefRcients in dS2 

We will begin with the (1+1) dimensional case for which the modes are 



(36) 
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So, the Fourier transform is: 

f{u) = I dt e^'^'e'^'^^" 



(,-ik/H 


POO 




'— oo 


^-ik/H 




V2k 





iiy/H 

r(-|)e--/2^. (37) 



Similarly, 



so that the modulus square of the coefficients of mixing in Eq. p4p are given by : 

That is the power spectrum per logarithmic band at negative frequencies (given by |/?iyp) is Planckian 
at temperature, H/2ti. At first sight this might look like the familiar result, well known in literature. 
However, there are some peculiar features which need to be commented on. 

Note that we have started with a solution to the wave equation in de Sitter background (given 
by Eq. ()36p ) and expanded it using the complete set of functions exp(ibii^t). These functions have 
no "legality" in the de Sitter spacetime since they are not the solutions of the wave equation. We 
could have, for example, used any other complete set of orthonormal functions in place of exp(ibiz^t) 
and could have defined the mixing coefficients through an equation like Eq. (j32p . The two properties 
which favour our choice are: (a) they were precisely the mode functions used to define the co-moving 
vacuum and (b) they are instantaneous, monochromatic, plane waves with respect to the co-moving 
time t. It is therefore interesting that the overlap between positive and negative frequencies in such 
an expansion gives rise to the thermal spectrum. 

If we had used the conformal time instead of co-moving time, then the result would have been 
very different — and very trivial. In terms of the conformal time, the modes are just exp(±zA:r/) at all 
77 and there is no mixing of the positive and negative frequencies defined with respect to ij. So if we 
had defined another set of mixing coefficients with an equation like Eq. (|32p but with conformal time 
r], then we would have got the trivial result = 0. So if we define the vacuum state with respect to 
conformal time and work entirely in terms of conformal time we will see no trace of thermal spectrum 
in the de Sitter universe. This is, of course, obvious from the fact that the metric is conformally 
flat in {r],x) coordinates and the scalar field theory is conformally invariant in (1+1) dimension; so 
we are back to the evolution of inertial vacuum in fiat spacetime. On the other hand, the modes 
undergo exponential redshift when frequencies are defined with respect to co-moving time and — 
as we had already mentioned — the exponentially redshifted wave will lead to a thermal mixing 
coefficient. We will next see that the situation is somewhat different in the (1+3) dimensional case. 

4.2 Mixing coefficients in dS4 

In this case, which we want to study in detail, it is convenient to work with a general mode function, 
having arbitrary coefficients Ak and B^- From Eq. (|26p . we have, 

Mt) = AkSkit) + Bkslit) = Vi'^(t) + (0 (40) 
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where 



Sk{t) 



exp 



ik 
H 



(1 



-Ht\ 



iH 

T 



+ e 



(41) 



Recall that taking Ak = 1 and Bk = gives Bunch-Davies state and for the co-moving state 
the corresponding values are provided by Eq. (j3ip . We again choose the normalization such that 
/k = Sk exp(zk • x) satisfies the standard orthonormality conditions (/k, /k') = (5^(k — k') with respect 
to the Klein-Gordon inner product: 



($i,$2 



(42) 



where the integral is evaluated over a hypersurface of constant t. 

Some amount of algebra yields the following results for the Fourier transform of the respective 
parts, V'fc^^(*) and ip^^\t): 



(2fc)3/2 



R 



r 



H 



iv 
H 



1 + 



H 



H 



(43) 



Taking the square of the modulus of the above expressions we get, 



2 f^e^" 










2 /3 






e/3^-1 





(44) 



where /? = 2ti/H. For negative frequencies, the forms of power spectrum are: 

,,2 



H/'"(-.)i^=i;iAf^ 



1 + 



1 + 



(45) 



Let us first consider the Bunch-Davies state, with Ak = 1 and = 0. This gives ai, = f^^\u) and 
(3y = f^^\—v) SO that we find: 



" 2k^u - 1 
1/3 ,2 = ^^ 



(46) 
(47) 



In contrast to the (H-l) dimensional case, these are not thermal, due to the extra factor + / H'^). 
So the expansion of the universe leads to a mixing of positive and negative frequencies but the 
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resulting mixing coefficients do not have a thermal form. It should be noted, however, that the ratio 
of the mixing coefficients is 

K\' |/W(z.)|2 ^^^^ 

When the field ipkit) couples linearly to a detector, the rate of upward and downward transitions 
between any two levels of the detector will be determined by the mixing coefficients. Therefore, 
when the condition in Eq. ()48p holds, one is led to a level population in the detector at thermal 
equilibrium with the temperature /3~^ = H/2-k. Any multiplicative function h{v'^) with and j3y 
drops off in the ratio. [Usually, one works with Bogoliubov coefficients which satisfy the constraint 
|ap — |/3p = 1; in that case, if Eq. ()48p holds, then |/3p must be thermal. In the case of mixing 
coefficients we have defined, the condition |ap — = 1 does not hold which allows extra factors 
like (l + iy^/H^).] 

Let us next consider the co- moving vacuum which holds more surprises. We now require the 
square of the modulus of complete f{—v) which is combination of individual quantities, \f^'^\—^)\'^, 
|/(^)(— z/)p evaluated above and a cross-term given by. 



2z.|/«(-i.)||/(2)(-z.)|cos( 



IP 



l3u/2 



■ COS 9 



(49) 



where 



9 = avg{f'>,f 



(1) f{2)^ 



(50) 



The complete expression becomes: 



1 + 



if2 



\Akf + IBkl'^e'^n N + 2\ Ak\\l3kWN{N + I) cos 9 



(51) 



where 



N 



1 



1 



(52) 



is the Planckian factor and Ak,Bk given by Eq. ()3ip . This result shows that for a massless scalar 
field prepared in the co-moving vacuum state, we obtain an expression having the Planckian factor 
with the temperature H/2tt. In addition, we obtain an interference term involving ^JN{N + V) 
which can be thought of as the fluctuation in the occupation number in thermal equilibrium. This 
factor has been noticed earlier [TT] in the case of horizon thermodynamics though no clear physical 
explanation is available. As far as we know, this has not been noticed earlier in the case of de Sitter 
spacetime in any context. 

5 Detector response in de Sitter spacetime 



The mixing coefficients defined through Eq. ()32p are directly related to the response of a co-moving 
geodesic detector in Friedmann universe. Since the clock carried by such a detector will measure the 
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co-moving time t, the rate of transition between the levels of the detector will involve the factors 
exp{zizitAE) where AE is the energy difference between the two levels. This gives an operational 
meaning to the mixing coefficients and we will show that the response of a co-moving, geodesic 
detector shows features very similar to what we obtained in the last section. 

Consider a stationary detector, located at the spatial origin in a de Sitter spacetime and coupled 
to the massless scalar field by monopole interaction. The amplitude for excitation of this detector 
during the time interval (— T, +r) due to its interaction with the scalar field can be computed, in 
first order perturbation theory as: 

jz/^^ = M j dTe'^^(lk|$(x[r])|0) (53) 

where M = iA(£'|m(0)|£'o) is amplitude of transition in the internal levels of the detector with A as 
the coupling constant and m(0) is the detector's monopole operator. (In the above expression, we 
are confining our attention to a final field state containing a particle with a specified momentum k. 
The total excitation probability for the detector is obtained by integrating |McP over all k.) The 
detector interacts with the field only during the period —T to T and x"(r) = x"(t) = (i, 0,0,0) is 
the trajectory of the detector. Expanding <I>(x[t]) as in Eq. ()14p . we find that the only term that 
survives in the T — t- oo limit is the negative frequency term. The amplitude arising from this term 
is given by, 

fT 



=24 = -A^ 



dt e 



iut 



Jk/H 



iH 

T 



+ e 



-ik/He- 



+ 



-ik/H (iJi ^ 



Ht 



Jk/He- 



This can be recast as 



where 



(54) 



(55) 



^^^i^t^-iki,/He- 



(poo r°° \ 

/ dt- dt- dt)e"''e-'^f'/'''~"\ 
J -00 J -00 Jt ) 



The above integral can be evaluated to give. 



g-W2^^gf InM 



where r(a, 6) is an incomplete gamma function. With this the amplitude becomes 



k 



At 1 



1 



(56) 



(57) 



(5^ 
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where we have ignored terms coming from differentiating the gamma functions since those are purely 
oscihatory and can be made to vanish in the large v limit by using the standard ie prescription. The 
probability Pk for the transition is now given by, 

= l^kl' =^ (i + ^) (lAI ViHP + l^.l 

^2\Au\\Bu\\h{y)\\n[-v)\ cos9) (59) 

To avoid the transients arising due to finite T, we will take the limit of HT S> 1. In this case, an 
elementary computation gives: 



nu/2H / ^"/H 



SO that 



\h{v)\ = ; \ cos 61 +0(6 ) 



u k \ V 

where N = {e^'^ — 1)~^ and /3 = 2it/H. Therefore, when HT ^ 1, we get the transition probability 
to be, 

= ^^^^ i^ + w)^ [(''^'^'' + l^fcl'e^"^) N + 2\Ak\\BuWN{N + l)co^e\ (62) 

A comparison with Eq. ()5ip shows that the detector response is triggered by essentially \f{—v)\'^ 
which should be obvious from the fact that the amplitude in Eq. (j53p picks out the negative frequency 
component of the field when the time integration is extended over the range (—00,00). This result 
shows that our mixing coefficients have a direct connection with the operational definition of particle 
content, as determined by the detector response. 

The above result is general and is valid for arbitrary Ak,Bk- By taking specific values we can 
determine the detector response in Bunch-Davies and co-moving vacuum. In the Bunch-Davies case, 
we have Ak = 1 and Bk = giving 

-^^(-^)^ 

This result shows that the detector response does pick up the extra factor (1 + v'^/H'^) just as the 
mixing coefficients do. (The same factor has been noticed earlier in ref. [12j). The corresponding 
result for co-moving vacuum can be obtained by substituting Eq. (j3ip into Eq. (j62p but the result has 
no special features worth mentioning. The above results arise because, by definition, the geodesic 
detector measures the co-moving time t. 
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6 Relation to the results in static coordinate system 



In Sec. [21 we briefly described how thermal nature of the de Sitter horizon arises in the static 
coordinate system and in the last few sections we studied the field theory in Priedmann coordinate 
system. Since both coordinate systems coexist in part of the de Sitter manifold, one can make an 
explicit comparison of the quantum states defined in these two coordinate systems. (This is similar to 
comparing the states in inertial coordinate system and Rindler coordinate system in fiat spacetime.) 
For this, we need to compute the relevant Bogoulibov coefficients on a spacelike hypersurface between 
the relevant mode functions by using the Klein-Gordon inner product. As we shall see, this is fairly 
straightforward in (1+1) but somewhat complicated in (1+3). 

6.1 Comparison in dS2 

We begin by noting that the metric, 

ds^ = dt' - e'^^'dx' (64) 
in (t, x) coordinates can be written in the static coordinates {t, x) as 

ds^ = (1 - H^x^) dP - (1 - H^i^y^ dx^ 

= (1 - H^x\x^)) {dP - dxi) (65) 

where 

x = e^^x] i=t-l/2H\n{l- H^x^) (66) 

and 

is the tortoise coordinate. Using these transformations we can express the field modes in the time- 
dependent dS2 coordinates in terms of the static coordinates. We will focus on a fixed {k > 0) mode 
so that the mode function 

fk{t,x) = _^e-'^l^e'^''e'^l^^~"' (68) 
v2/c 

becomes 

fk{t,x) = ^=exp -J- H -J- = ——e^y^'^i'' (69) 

in static coordinates where u = t — x^. In the static deSitter patch, conformal flatness of the metric 
in Eq. ()65p allows us to write down the solution to the field equation as exp(±icjM). This allows the 
expansion 

= ^ (ftc.e^^'^" + 6ie^'^") (70) 



/2a; 

etc. which is valid on the complete manifold. We now need to determine the Bogoliubov coefficients 
that relate the above two sets of field modes. These are given by the standard Klein-Gordon inner 
product. 



\j_d5:.{^^difk-fkdi^^) (71) 
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where the integral is over any spacelike hypersurface. Choosing t = surface, the above integral 
over x<t: can be recast as: 

f3^^ = -= du [e-^'^^dufk - fk due-'-'') (72) 
Integrating the first term by parts gives, 

/oo 
dne"^""/fc(^^) + /fce 
-oo 



iOIMlOO 

I ~oo 

OO 
OO 



= V2uj I dne"*'^"/fc(u) (73) 
since the second term vanishes. Thus 



oo 



We find that modulus |/3i^fc| is again Planckian at temperature, H/2'k: 



This shows that the Bunch-Davies vacuum (which is the same as the co-moving vacuum in (1+1) 
dimension) has a thermal character in the static patch bounded by the horizon. 

6.2 Comparison in dS4 

The transformation from the Friedmann coordinates to static coordinates goes through in dS4 exactly 
in the same way as dS2. The metric 

ds^ = dt^ - e^^* {dr"^ + r^dn^) (76) 

in (t, r, O) system can be written in the static coordinates (t, f, Q) as 

ds^ = (1 - H^r^) dP - (1 - H^rY' dr^ - f^dQ^ 

= (1 - H'^r\r,)) {dp - drl) - r\n)dn^ (77) 

with the same transformations as before 

r = e-^V ; i=t- 1/2H In (l - if^f^) (78) 

and defining the tortoise coordinate 



df 

(1 - H^f^) 



(79) 
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In the static coordinates, the field equation reads 



$(t,f,O) = (80) 



where /(f) = (1 — H'^r'^) and L is the standard angular Laplacian operator. Taking 
$ = (j)i{f)Yim{^)e~^^^ / f, we find that (pi{r) satisfies the equation: 



- oj'c^i - ^-p r'f- {^))- ' 1^ 0^ = 0. (81) 
r dr \ ar \ r J J 

Since /(f) vanishes at the horizon f = 1/H, only the s-mode makes a dominant contribution near 
the horizon and hence we will focus on / = mode. For this mode, the wave equation becomes 

where the prime denotes derivative with respect to f. Clearly, in the near horizon limit (/ — )• 0), the 
solutions behave as exp(ibia;r,,). Thus near the past horizon, f — )• 1/H and t — ?• — oo, the modes in 
the static coordinate system behave as exp(zbia;t)). 

On the other hand, the modes describing the Bunch-Davies vacuum can be expressed in spherical 
coordinates by the standard plane wave expansion 



^^"^ = ^ E^'(2^ + m{km{^ose)e^''^^-"' + e-^*j (83) 



Using the transformations in Eq. (j78p we can express this in {t, r) coordinates. Concentrating on 
the s-wave contribution we obtain 



/2fc V / V ^ 

which, near the past horizon, f — t- 1/H and t — t- — oo, behaves as 



2k 

We now use the fact that the Klein-Gordon inner product between the field modes is independent 
of the surface over which it is evaluated. It is, therefore, convenient to evaluate the Bogoliubov 
coefficients on a space-like surface very close to the horizon. Since the Bunch-Davies mode behaves as 
^t{k/H)e " while the static modes behave as exp(±ia;u), it is obvious that the Bogoliubov coefficients 
defined in Eq. ([75]) will give 

/"OO 

fi^k = y/2^ dve-'^^^f''\v) (86) 



which has a thermal character: 



,2 P a 2^ 
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We again see that the Bunch-Davies vacuum has a thermal property when viewed in the static patch 
in (1+3) dimensions as well. In this sense, (1+1) and (1+3) dimensions behave identically. It is 
also straightforward to show that a detector at rest in the static coordinates will perceive a thermal 
radiation in the Bunch-Davies vacuum state. On the other hand, we saw earlier that a freely-falling 
detector will also see the modified thermal spectrum [see Eq. ()63p ] in the same vacuum state. It 
should be noted that this is somewhat contrary to the results in black hole spacetime. 

Finally, we quote the result for the co-moving vacuum transformed to static coordinates. The 
analysis is again straightforward when we use the fact that the co-moving modes can be expressed 
in terms of the Bunch-Davies modes by the relation 

^l.f''\t)=A,i.^^^\t) + B,^(^^r{t) (88) 

Therefore, 

^iCM) ^ 1 U^^k/He-«^' ^ Q^^-^me~^^\ . (89) 

v2A; ^ ' 
on the past horizon. It follows that the spectrum is now given by 

k\^u.k? = (lAl' + l^fcl'e'^") /3iV + |A||Sfe|/3V^^(iV + l)cos0 (90) 

where 



We once again see that the co-moving vacuum introduces an interference term in the form of 
\JN{N + 1) even when we compare the modes between Friedmann description and static description 
suggesting that there must be some physical explanation for the origin of this factor. We hope to 
address this question in a future publication. 



7 Quantum fields in quasi-de Sitter spacetime 

So far, we have been concentrating on the features which are special to de Sitter spacetime. However, 
in the evolution of the real universe, it is impossible to obtain a pure de Sitter evolution due to the 
presence of external matter. Both, during the inflationary phase as well as during the late time 
acceleration phase, we only have a quasi-de Sitter phase rather than a pure de Sitter universe. 
In this section we will extend the formalism described earlier to a quasi-de Sitter spacetime by 
determining an approximate solution to the wave equation. This approach is quite general and can 
take into account any first order deviation from the pure de Sitter universe. After developing the 
formalism we will apply it to a specific example to illustrate its utility. 

7.1 The perturbative framework 

Consider a Friedmann spacetime with the scale factor given by: 

a(t) = e(-f^*+^^(*)) « e^' (1 + eX{t)) = ao + eXao (91) 

which can be treated as quasi-de Sitter if the condition A <^ XH is satisfied. In the above expansion, 
we have retained the perturbation to first order as indicated by the bookkeeping parameter e (which 
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will be set to unity at the end of the computation). Correspondingly, the mode functions, which are 
the solutions to the wave equation in the perturbed metric, will differ from those in the de Sitter 
spacetime by a small amount: 

il^{t) = Mt) + e6^{t) (92) 

where ipoit) is the unperturbed mode function and we have omitted the subscript k for notational 
simplicity. Substituting the above expressions for a{t) and ijj{t) into the time dependent part of the 
wave equation written in the form: 

+ 3 - -^ + -:7V' = 0. (93) 



dt"^ \a ) dt a? 

we get. 



— 2— + 3i^— — + 3A-— - 2^A'0o + ^'JV' = 0. (94) 
at^ dt dt 



This equation can be solved by writing 5ip = ipos. The function s{t) then satisfies the equation. 



fi{t) = 2—X-3X-f- (96) 



where 

!^A-3A^ 

acts like a source term. Eq. (j95p is first order in ds/dt and hence can be immediately integrated. 
(This result holds for a generic class of second-order homogeneous linear differential equation; see 
Appendix [B] for details). The solution for s{t) is 

s{t) = C J' dtVo~'e-'^*' + dt' Vo"'(t') j' di'Vol^'Oe'^'^'Vli") (97) 

where C is a constant of integration. Thus, given a model for A(t) and appropriate boundary 
conditions we can, in principle, solve for the perturbation Sip by this method. 

7.2 An example: late-time accelerated phase of the universe 

As an illustration of the above method, let us consider the late time accelerated phase of the universe 
containing dust-like matter and a cosmological constant. The expansion factor of such a universe is 
given by 

/ o \ 2/3 

a{t) = 2^/3 sinh -Ht (98) 



V 2 

In the spirit of the above discussion, we will treat this as a perturbation to an exact de Sitter universe 
and write 

a{t) « e^* (l - ^6-3^*] (99) 
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where A(t) = —{2/3) exp{—3Ht) which vanishes as t goes to infinity. This behaviour suggests that 
we use the boundary conditions s(oo) = and s(oo) = in our general sokition given by Eq. ([9] 
In the pure de Sitter case, the mode functions can be taken to be 



1 

v2k 



^ + e"^* ) (100) 



which amounts to taking = 1 and i3k = in the Eq. (j26p i.e., we have chosen to work with 
Bunch-Davies state. Calculating the integrals in Eq. ()97p is straightforward (see Appendix [C] for 
details) and we obtain, 

sit) = r^e-^^^ (101) 
Therefore, the first-order change in the mode function is given by: 

6m = Mt)s{t) = '-^^ (l - §e"^*) e"^^*et-"\ (102) 

We can now compute the Fourier transform of this expression to determine the first order correction 
in Fourier space: 



Sf{u) = / dt6i;{t)e 



tut 



lbHy/2k J-oo V H 

7ik 



/OO 



W2H(g_^ r 5-^) (103) 



The resulting power spectrum, to the lowest order, is given by: 

z.|F(z.)|2 « z.|/(z.)|2 + ei.Re[2riu)5f{i^)] (104) 
where we have reintroduced e for bookkeeping and 

.Re[2r (.)5/(.)] = Jt^T (i44^ + 64^ - 79^ + ^ j (105) 

is the correction to the power spectrum in the case of quasi-de Sitter phase arising from the matter 
contribution in the late-time acceleration. 

8 Conclusions 

The periodicity in the Euclidean time allows us to attribute a temperature H/2tt using the static 
coordinates on the de Sitter manifold. In this sense, de Sitter spacetime behaves just like other static 
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spacetimes with a horizon. However, in such an analysis, one has to define a vacuum state on a 
T = hypersurface in the Kruskal-hke coordinate system which is not static. In this particular case, 
thermal nature of the de Sitter horizon arises because the vacuum state in Kruskal-type coordinate 
system leads to a thermal density matrix for the observers bounded by the de Sitter horizons. 

The de Sitter spacetime is unique in the sense that it also allows introducing Friedmann coordi- 
nates in which the metric is homogeneous. This, in turn, reduces the field theoretic problem to that 
of a quantum oscillator with a time dependent frequency. In such a, time dependent, background 
there is no unique definition for the vacuum state and the best one could do is to introduce well 
motivated vacuum states and study their physical properties. Quite generically, such states can be 
introduced by giving a suitable boundary condition for the mode functions at some time t = to. The 
question arises as to whether one can understand the thermality of de Sitter universe working en- 
tirely within the Friedmann coordinates i.e., without comparing the results between Friedmann and 
static coordinates. (We have not seen such a derivation in the literature for a massless scalar field.) 
We investigated several aspects of this question both in (1+1) dimension and in (1+3) dimensions 
in this paper. 

Two natural vacuum states one can introduce are the Bunch-Davies and co-moving vacuum states 
in this spacetime. In (1+1) dimension, the Bunch-Davies vacuum state corresponds to choosing the 
modes to be positive frequency with respect to the conformal time r] in the asymptotic past while the 
co-moving vacuum state corresponds to imposing the positive frequency condition at some arbitrary 
instant of time t = to- It turns out that both these states are identical in (1+1) dimension. To study 
the time evolution of this state, we expand the mode function in terms of positive and negative 
frequency modes defined with respect to the co-moving time. The mixing of positive and negative 
frequency modes then reveals a thermal character with temperature H/2tt. This, of course, does 
not happen during the time evolution in the conformal time; the positive frequency mode remains a 
positive frequency mode at all times. 

The situation in (1+3) dimensions is quite different. To begin with, co-moving vacuum state and 
the Bunch-Davies vacuum state do not coincide in (1+3) dimensions. Further, the mixing coefficient 
between positive and negative frequency modes does not have a pure thermal character (and is 
modified by an extra frequency dependent factor) in the case of Bunch-Davies vacuum. The result 
for the case of co-moving vacuum is more complicated and involves an interference term containing 
y^/V(iV + l) factor which is reminiscent of the fiuctuations in the occupation numbers of massless 
thermal radiation. 

The physical meaning of the mixing coefficients introduced to analyse the above phenomena 
can be understood by studying the response of particle detectors in the de Sitter spacetime. We 
computed the rate of excitation of a geodesic detector evolving in co-moving time. This rate exactly 
matches with the particle content of the state as determined by the mixing coefficients in both 
Bunch-Davies vacuum and co-moving vacuum. 

We also compared the states defined using Friedmann coordinate system with those defined 
using the static coordinate system. This requires evaluating the necessary Bogoliubov coefficient 
between the mode functions defined in the static patch and Friedmann patch in the region of the 
manifold where they coexist. We found that the Bunch-Davies vacuum appears to be a thermal 
state for static observers bounded by the horizon, both in (1+1) and (1+3) dimensions. This is in 
contrast with the results obtained within the Friedmann coordinate system where the results for 
(1+3) dimensions differ from the results for (1+1). On the other hand, the co-moving vacuum in 
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(1+3) dimension, defined in Friedmann coordinates, does not have a simple thermal interpretation 
in the static coordinates. 

In the last part of the paper, we studied the effects of small deviations from de Sitter evolution 
and the resulting corrections to the mode functions. This formalism is sufficiently general to handle 
any functional form of the deviation in the lowest order of perturbation theory. As an illustration 
of this formalism, we studied the deviations in the power spectrum arising due to the existence of 
pressure-free matter during the late time accelerated phase of the universe. This formalism might 
have applications for studying the spectral deviations in the case of inffationary universe as well. 
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A Calculation of the Fourier transform in Eq. (1331) 



To evaluate 



we define u = e and b = kfi/H. This gives, 

"OO 



H 



duu He 



ibu 







^exp 



H 



H\H 



(106) 



B A result in perturbation theory 

Consider a generic second-order homogeneous linear differential equation 

a{t)x{t) + h{t)x{t) + c{t)x{t) = 

Let xq be the solution of above equation for some functions ao(i)) ^o(*) and co(t). We are now 
interested in the corresponding solution of the equation when the parameter functions a, h and c are 
perturbed about their original forms. Then, to the first order in perturbation, we have 

a^Sx + x^Sa + h^Sx + xo6b + cqSx + xq6c = 0. (107) 

Scaling the perturbation 5x with the unperturbed solution as, 5x = xqs, gives for s{t) the equation 

s + A{t)s = B{t) (108) 
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with 



A{t) 
B{t) 



xo ao 

xq 5a xq 6b 6c 
xq ao Xq ao ao 

which is a first order differential equation in s and can be solved immediately. 



(109) 



C Solution of Eq. (1971) for late-time accelerated phase 



The basic ingredients that go in are 

1 







exp 



ik 
H 



l-e 



-Ht\ 



iH 



+ e 



-Ht 



which gives 



and 



2^-2Ht 



■00 — ik'^e 
^ ~ H{i + |e-^*) 



fi{t) = 2k^e~'^"^\ - 3A 



;0o 



4 , „5^, Gik^e-'^' 
'3 ^{i + k/He-Ht) 



so that under the conditions s(oo) = s(oo) = 0, we can set C, the constant of integration in the 
homogeneous part to be zero and obtain: 



e H 



-3Ht' roo 



, // 2ik ^ — Ht ou+'f 
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44 _ 4 — e 



~2Ht'' 



ik 

+ 10-e 



-Ht" 



k 

57 ( 3, -2i-^e 



-3Ht' 



-147 (2,-2z|:e-^*')+7(4,-2zAe-^*' 



-Ht' 



Noting that, 

-147 (2, x)+-f (4, x) - 57 (3, x) = -18 + e^(18 - 18x + 2^^ + x^) 
we can evaluate s(t) in the late-time approximation to give 



sit) 



-5Ht 



(110) 



as the leading order correction term. 
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